Emission Line Flickering from the Secondary Star in Cataclysmic 
Variables? A study of V3885 Sagitarii 
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ABSTRACT 

Spectrophotometric observations of Ha and He I 6678 emission lines of the nova-like Cata- 
clysmic Variable V3885 Sgr are presented and analyzed. The binary orbital period was determined 
as P = 0.20716071(22) days. Doppler Tomography was performed with both Ha and He I lines. 
Disc radial emissivity profiles were also computed. The tomography mapping of flickering sources 
was performed using the Ha line, from which we concluded that the flickering is not uniformly 
distributed on the disc. The observed tomogram of the flickering was compared with simulations, 
suggesting that the most intense flickering source in the Ha is not located in the accretion disc. It 
is proposed that the main line flickering source may be associated with the illuminated secondary 
star. 



Subject headings: accretion, accretion disks - 
(V3885 Sgr) 

Introduction 



binaries: close — cataclysmic variables — stars: individual 



Cataclysmic Variables (CVs) are close binary 
systems composed of a white dwarf, also called 
primary, and a red dwarf star, the secondary. The 
secondary fills its Roche lobe and transfers mat- 
ter to the primary. Due to the binary angular 
momentum, the matter does not fall toward the 
white dwarf, forming an accretion disc. The accre- 
tion disc emission is complex, it has regions of en- 
hanced emission like the hot spot, the locus where 
the stream of matter coming from the secondary 
reaches the disc, as well as the boundary layer, the 
inner part of the accretion disc where the matter is 
decelerated to the white dwarf equatorial velocity. 
The disc's vertical profile is not flat; it presents a 
higher thickness in it s outer parts. A review about 
CVs is presented by Warner ( 19951 ). 

The CVs are grouped into many classes and 
subclasses; one of these is the nova-like class 
that harbor non-eruptive objects. V3885 Sgr is 
one of the brightest CVs. It is a non echpsing 
system classified as an UX UMa nova-like star. 



V3885 Sgr was discovered by iBidelman &: Bond 
(Il968r) w ho detect ed br oad Balmer emission lines. 
IBond fc L andolt ( 19711 ) identified the continuum 



slope as simi lar to those of DA white dwarfs. 
Wegner (|l972h observed emission lines superposed 



on absorptio n lines identifying the object as a 
nova-like CV. Hesser et al. I ((l972l) first noticed the 
occurre nce of rapid phot ometric variability in the 
system. ICuinan fc Sim] (|T982) carried a spectro- 
scopic study in the UV noticing the existence of P- 
Cygni profile in the C IV line, which indicates the 
presence of a wind in the system. Fr om UV , opti- 
cal and IR continuum distributions, Haud ( 1987 ) 
described the spectral energy distribution of the 
disc as one o f an stationary disc with high mass 
transfer rate. iHartlev et al. ( 2005 ) calculated Por b 
— 0.207135(15) d for this system, solving previ- 
ous inconsistencies concerning the orbital period. 
They also presented a radial velocity study, H/?, 
H7 and Hel 4471 Doppler maps, and they showed 
evidences of irradiation for the secondary star. 

Flickering or rapid variability is observed as 
a stochastic variation of the emitted radiation, 
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with time-scales from seconds to tens of min- 
utes and amplitude ranging from a few percent 
of magnitude to more than one magnitude. Flick- 
ering is associated with mass transfer and oc- 
curs in objects other than CV s, like symbiotic 
stars ([Mikolaie wski et alj 
naries (je.g 



1990 



and X-ray bi- 
Furthermore, 



■ g. Baptista et all 12002 
flickering does not only occur in systems where 
an accretion disc is present but it is also observed 
in magnetic CVs. The flickering in the contin- 
uum cou ld be produced in differen t regions of the 



system. IWarner fc Natherl (|l97lh observed that 



the flickering in U Gem disappeared during the 
eclipse, and since the central part of the disc is 
not occulted during the eclipse, they associated 
the fli ckering source with the hot spot. IVogt et aT 



(|1981I ) verified that the flickering in OY Car per- 
sists during the hot spot eclipse. The sani e was 



PattersonI (Il98lh and 



observed for HT Cas bv _ 

Z Cha by IWood et all ()l986l ). iHorne fc Stienine 
identified the flickering in some eclipsing 



systems as being originated n ear the cente r of the 



accretion disc. In addition, Horne et al. (1993) 



associated the flares obser ved in OY Ca r with the 
innermost part of the disc. Bruch ( 1996h observed 
that the flickering forming region in Z Cha is lo- 
cated near the white dwarf and that there could 
be other flickering sources seen in different photo- 
metric states. 

In the next section we detail the observations 
and data reduction procedure. In section 3 the 
average spectrum is described and a long term 
ephemeris is given. The system masses and or- 
bital inclination are constrained. Doppler Tomog- 
raphy and Flickering Doppler Tomography results 
are presented. Furthermore, the observed flicker- 
ing tomograms are compared with flickering sim- 
ulations. The implications of our results are dis- 
cussed in section 4. The conclusions are outlined 
in section 5. 

2. Observations and Data Reduction 

The time-resolved spectroscopic observations 
were made using the 1.60-m telescope at Ob- 
servatorio Pico dos Dias (LNA/CNPq) and the 
1.5-m telescope at Cerro Tololo Inter- American 
Observatory. For more details about the obser- 
vations, see the observations journal (Table 1). 
The atmospheric conditions on the observing runs 



ranged from photometric to spectroscopic. Dif- 
ferential spectrophotometric techniques were em- 
ployed. The LNA observations were made using 
a Cassegrain spectrograph, equipped with a 1200 
grooves/mm grating blazed at 4550 A, and a back- 
illuminated 1024x1024 CCD. The CTIO observa- 
tions were also carried out with a Cassegrain spec- 
trograph, a 1200 grooves/mm grating, but blazed 
at 7500 A, and a 1200x800 Loral IK CCD. All 
the observations were centred at the Ha line. A 
spectral resolution of about 2 A was achieved. The 
spectral coverage also included the He I 6678 spec- 
tral line. 

Bias and flatfield images were taken to adjust 
the zero level count and to correct for response. 
When the observational conditions allowed, twi- 
light flat images were also taken for deriving the 
CCD illumination correction. The data were bias 
subtracted, flatfielded and illumination corrected 
using standard IRAlIl] procedures. The spectra 
were e xtracted usin g an optimal extraction algo- 
rithm ( Hornd 19861 ). for both the object and a 
comparison star. Each spectrum was then cor- 
rected for possible cosmic rays by interpolating 
neighbor pixels. We set the optimal extraction 
to trim only the most intense cosmic rays, as a 
sharp line profile feature (i.e. due to flickering) 
could be confused with a cosmic ray and trimmed. 
So we prefer to eliminate the less intense cosmic 
rays by hand. The cosmic rays in the sky back- 
ground regions were all automatically removed. 
The wavelength calibration was performed by in- 
tercalating the Ne-Ar arc lamp spectra with obser- 
vations of the variable star, resulting in an accu- 
racy of ~ 0.05 A RMS or 2 km s"^ The interval 
between arc-lamp exposures was chosen to min- 
imize the spectrograph mechanic flexure effects 
over the wavelength sca le. Wide slit tertiary stan- 
dard photometric stars (jHamuv et a l. 1994) obser- 
vations were taken for flux calibration. To allow 
the correction of light losses at the spectrograph's 
slit and to perform small sky transparency cor- 
rections, a nearby comparison star was included 
in the slit. Wide slit (7 arcsec) spectra of the 
comparison star were also acquired. This differen- 
tial calibration procedure also corrects the spectra 



^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation 
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from atmospheric dispersion effects. 

As the Ha hne region includes telluric lines and 
bands, we have performed a correction using a 
telluric template with the same spectral resolu- 
tion. The telluric template was scaled and shifted 
to correct the telluric features using an averaged 
spectrum of the target. This correction was then 
applied to each spectrum of our dataset. For this 
reason, the correction should be good enough to 
the average spectrum and may not be so efficient 
for each individual spectrum. Individual correc- 
tion factors for of each spectra are unreliable be- 
cause of their low S/N ratio. All derived quantities 
involve averaging data on phase bins therefore we 
do not foresee any impact on the measured quan- 
tities. 

3. Results 

3.1. Spectral Features 

The averaged spectrum (fig. 1) has a spectral 
coverage of 6400-6714 A, showing Ha and He I 
6678 emission lines. The continuum has a blue 
slope, without any other visible feature. Ha ap- 
pears as single-peaked and He I as double-peaked. 
The former line presents a small absorption com- 
ponent at the blue side of the profile, indicat- 
ing a P-Cygni profile or an underlying absorp- 
tion. The Ha line has a FWHM of 780 km s'^ , 
EW ~ —4.1 A, and the line profile integrated flux 
F is 5.7 10^13 erg s^^ cm^^ ^^lile the He I 6678 
line has a FWHM of 820 km s'^, EW 
and integrated flux F = 6.4010"^'* 
thus the equivalent width of Ha line is 8.5 times 
the He I one. 

It was not possible to identify any absorption 
spectral feature from the secondary in the phased 
spectra or in the averaged spectra. So, the sec- 
ondary radial velocity amplitude if 2 and the com- 
panion spectral type could not be directly deter- 
mined in the present work. 



erg s ^ cm ^ 



3.2. System Parameters 



A long term ephcmeris was derived using un- 
ambiguous cycle-counting between all observation 
nights in our dataset. The calculated orbital pe- 
riod for V3885 Sgr is Forb = 0.20716071(22) d and 
epoch Eq = 51942.08371(46) MHJD, defined as 
the positive to negative crossing of radial velocity 



curves for the line wings. The uncertainties were 
obtained from a least squares linear fit to all cy- 
cle timings. The residuals in the 0-C diagram are 
well behaved and there is no suggestion of a pe- 
riod derivative. This period is in a greement with 



the 0. 207135(15) d period found bv lHartlev et al 
(|2005h within a 2a interval 



Using the lSchneider fc Youngl ( 1980[ l method to 
measure the velocity in the line wings, we have 
constructed a diagnostic diagram for both Ha and 
Hel as a function of the convolution mask Gaus- 
sians separation a. The primary radial velocity 
(Ki) values remained almost constant inside the 
interval 300 < a < 600 km s^^, presenting a 
slightly decrease in velocity for a > 600 km s~^. 
The relative error on Ki as a function of a showed 
a minimum at a « 500 km s^^. The RMS ex- 
hibited this same behavior. From this diagnostic 
diagram, looking for a range in a that presents 
a plateau in K and minimum relative error and 
RMS, we estimated the primary radial velocity 



from KHa = 139(25) 
sion line and KheI = 



km s for the Ha emis- 
188(25) km s'^ for Hel. 
These uncertainties were estimated on the basis of 
the assumption that the radial velocity measure- 
ments from both Ha and Hel should be compati- 
ble. The Ki value does not vary significantly with 
the Gaussian mask width. The Gaussian masks 
FWHM widths were 200 km s"^ for Ha and 150 
km s~^ for Hel. These radial velocity values are 
compatible wi t h K = 166(6) km s~^ presented by 



Q 4g ^ Hart lev et al. ( 20051 ). The systemic velocity was 

also derived for both lines, jHa = -43(25) km s~^ 
and ^Hei = -84(25) km s~^ 

The double-Gaussian convolution technique for 
radial velocity me asurements is not immun e to 
systematics errors (jUnda-Sanzana et al.ll2006r ). In 
some cases it produces wrong masses estimates, 
even with large datasets like the one present here. 
Unfortunately, there is no other way to constrain 
the primary mass in this non-eclipsing nova-like 
CV. 

The mass diagram for V3885 Sgr (fig. 2) was 
constructed using the Ha radial velocity semi- 
amplitude and its error. The continuous verti- 
cal line at the right indicates the white dwarf up- 
per mass limit. The two vertical lines at the left 
are the white dwarf lower mass limit for i = 75° 
(dashed line) and i = 45° (continuous line). These 



limits were obtained from each line FWZI, consid- 
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ering a Keplerian disk and a ma ss-radius relation- 
ship (|Hamada fc Salpeteii 119611 1 for the primary 
star. The diagonal straight line is the equal mass 
line, plotted only as a reference as the stable mass 
transfer line is not straight along the mass dia- 
gram and it is also model dependent. The quasi- 
horizontal continuous line represents the upper 
mass limit for a main sequence secondary filling 
its Roche lobe. The Mi,M2 relation is plotted for 
a few inclination values (i = 35° is plotted as the 
dotted lines, i = 45° as the dashed lines and i = 75° 
as the dot-dashed hues), considering Ki and its 
uncertainty. The region of most probable masses 
is shown in gray. The primary mass is uncon- 
strained between 0.3 Mq and approximately 1.4 
Mq, and the secondary mass is possibly between 
0.2 and 0.55 Mf7^. The prim ary mass range ob- 
tained bv lHartlev erall (|2005l ) is 0.55 < Mi < 0.8 
Mq. Another parameter that can be constrained 
from the mass diagram is the orbital inclination, 
which is found between 45° and 75°. The analysis 
and conclusions in the following sections are not 
dependent on a precise Mi value. 

3.3. Doppler Tomography 

The Doppler Tomog raphy method, proposed by 
Marsh fc Ho"m3 (|l988h . uses the velocity mforma- 
tion from spectral line profiles at various orbital 
phases in order to reconstruct the binary system's 
emissivity in the velocity space. The method is 
based on a few hypothesis listed below. All the 
emission line production regions must be always 
visible and their line optical depth must be neg- 
ligible, which is violated if observations during 
eclipses are used. The line flux at any point of 
the binary rotating frame of reference must be 
constant. Outbursts represent a flux variation, 
breaking this hypothesis. All the motions must 
be parallel to the orbital plane. Wind originating 
from the disc has a velocity vector out of the or- 
bital plane. All the velocity vectors rotate with 
the binary system. The last hypothesis is that the 
intrinsic line profile width in any point is negligi- 
ble; saturation effects at large optic al depths can 
widen the fine profiles (|Hornelll995l ). 

The Doppler tomograms from the Ha and He I 
line profiles were obt ained using the filtered back- 
projection method (jRosenfeld fc KakI 119821 ) and 
are shown in Fig. 3. The coordinate system is 
defined with the x-axis pointing from the primary 



to the secondary and the y-axis pointing in the 
direction of the secondary's motion. Considering 
a Keplerian disc, the outer part of the tomogram 
(higher velocities) represents the inner region in 
the position space, and vice-versa. The Ha and 
He I tomograms were calculated using 1888 inde- 
pendent spectra. A 70 km s^^ FWHM Gaussian 
convolution filter was applied to the Ha line and 
a 120 km s~^ one was applied to the He I line. 
The tomogram resolution is a combination of the 
spectral resolution, convolution filter width and 
phase sampling. As the V3885 Sgr masses and 
orbital inclination cannot be strictly determined, 
the marked positions on the tomograms are plot- 
ted just as a reference. Distinct values of orbital 
inclination and masses would produce a rescaling 
of the marked positions. 

In the Ha tomogram (Fig. 3a) one can see 
that the most intense emission is located at the 
{—Vx,+Vy) quadrant, and we could not affirm 
based only on this map whether this emission is 
originated from the illuminated secondary or from 
the hot spot. It is not aligned with the position of 
the stellar centres or clearly over the stream tra- 
jectory. The stellar positions could move along the 
y axis direction depending on the actual orbital in- 
clination. An extended region in the {—Vx,—Vy) 
and {+Vx, —Vy) quadrants can also be seen, corre- 
sponding to disc regions opposite to the secondary 
star. In order to show more clearly these tomo- 
gram anisotropics, a radially symmetric disc image 
was subtracted from the maps (Fig. 3). 

The He I tomogram (Fig. 3b) presents a ring 
shape when compared with Ha, suggesting that 
the emission from the former is mostly produced 
in the disc itself. The spot in the {~Vx, + Vy) quad- 
rant on the Hel map is approximately at the same 
location as the Ha map spot. The asymmetric 
component is more structured in He I than in Ha, 
and clearly shows the enhanced emission in the 
{+Vx,—Vy) quadrant. 

Comparing the Ha and He I tomograms (Fig. 
3), one can see that the He I emission is more in- 
tense in the {+Vx, +Vy) quadrant and in the lower 
part of the {—Vx,—Vy) quadrant. The origin of 
such behavior is still unclear. It can also be seen 
that the Ha emission is more intense than the He 
I one in the {—Vx,+Vy) quadrant. The emission 
from the {+Vx,~Vy) quadrant in both lines can 
be explained if the stream does not hit the disc 



4 



in the {—Vx,+Vy) quadrant, but it flies above the 
disc, reaching the disc in the (+14, —Vy) quadrant, 
before the point of minimum approach to the pri- 
mary. Therefore, part of the stream material is 
not incorporated to the disc in the hot spot re- 
gion, but follows a free particle trajectory above 
the disc until it reaches the {+Vx, —Vy) quadrant. 
This hypothesis is plausible if one reminds that 
the accretion disc geometry is not flat, and so the 
height in the disc's outer rim is greater than in 
the central regions and, in addition, the primary 
gravity has an increasing vertical component for 
smaller radii. 

The consistency of the Doppler maps was ver- 
ified by comparing trailed spectra generated from 
the projection of the maps at many orbital phases 
and the trailed spectra generated from the ob- 
served line profiles (Fig. 4). The comparison be- 
tween the phase maps suggests a reasonable con- 
sistence of the tomographic reconstructions. A sig- 
nificant difference between the maps is related to 
the sharp and low velocity feature seen in the ob- 
served phase maps. This feature appears once per 
orbit in the observations while it is spread into two 
features on each orbit in the projected map. As 
the tomography algorithm considers that all the 
structures are visible all the time, this sharp emis- 
sion should be related to the illuminated face of 
the secondary star (see discussion in section 3.5). 

Emissivity profiles in velocity space can be 
obtained directly from the Doppler tomograms. 
That was done using Doppler tomograms for Ha 
and He I centred at the primary's position (Fig. 
2b). The emissivity profiles are similar to line 
profiles but they are free from the projection of 
intrinsically high velocity emission that may ap- 
pear at lower velocity in the line profiles. In order 
to produce the velocity profiles, the intensity along 
each concentric ring centred at the primary were 
represented by its mode. With this procedure, we 
verified that the Ha tomogram is more homoge- 
neous than the He I one, with the mode, median 
and mean yielding similar results. 

Doppler tomograms are always plotted in the 
velocity space, but one can transform the maps 
to position space, assuming a Keplerian velocity 
regime. One problem that could occur in this 
transformation arises from the fact that the disc's 
most internal regions (seen in the line far wings) 
are brighter and carry the larger uncertainties in 



the inversion process. To study the disc emissiv- 
ity profiles, the emission at a given radius was es- 
timated, obtaining the disc radial emissivity pro- 
files for both Ha and He I (Fig. 3a). This value 
is calculated as the mode along each ring because 
it is less susceptible to local variations due to the 
disc emission anisotropies. To obtain the abso- 
lute radial emissivity profiles, we used a distance 
to the system equal to 11 0(23) pc from Hippar- 
chos parallax ( ESAI IQQtI) and a colo r excess of 
E{B-V) of 0.02 ()Bruch fc EngellflQOl . The pro- 
files are shown in Fig. 5. The radial emissivity 
profiles may be fitted by power laws, Iho. oi r~^-^ 
and IhcI ct r~^-^. In this system the He I emis- 
sion seems more concentrated in the inner disc 
regions than the Ha one. The obtained radial 
emissivity indexes are similar to t h ose o btained 
for V841 Oph by iDiaz fc Ribeirol (|2003h . Ra- 
dial emissivity profiles were also obtained from 
Dogpler maps fo r CU Vel (|Mennickent fc Diaz 
19961), V347 Pup (iDiaz fc Hubenvl ll999l). RR Pic 



(IDiaz fc Ribeirol 120031) and U Gem during out- 



burst (jMarsh fc Horneiri990f ) 



3.4. Flickering Tomography 

We have explored until this point the "station- 
ary" line emission from V3885 Sgr. The line emis- 
sion intrinsic variability can be conveniently stud- 
ied using the tomographic reconstruction of the 
fiickerin g emission. In the flickering tomography 
method (|Diaal200l[ ). we start calculating the vari- 
ance of the observed line profiles into phase bins. 
Then, the intrinsic line variability is isolated from 
the instrumental, orbital, and secular components. 
The instrumental component of the variance is 
calculated as a combination of the Poisson noise 
and a Gaussian readout noise. The orbital com- 
ponent contribution is considered negligible, since 
the data is sampled into small phase bins. The 
secular component is the remaining long term vari- 
ability of the system. The effect of such a compo- 
nent was found to be negligible by replacing the 
long term average profiles by profiles computed 
from individual runs. Doppler tomograms are cal- 
culated from these intrinsic variance line profiles 
mapping the line emission variability in the sys- 
tem. The flickering tomography method is very 
demanding on the amount of collected data and 
also sensitive to the time sampling. 

Pure photon and readout noise were removed 
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from the data using the empirical noise model de- 
scribed above. However, residual noise from other 
sources (such as continuum subtraction and dif- 
ferential calibration errors) are still present in the 
variance data. 

A total of 1863 spectra from our original data 
set were used in the flickering tomography study. 
Some spectra were not used in this analysis be- 
cause they deteriorate the statistics on the quan- 
tities calculated with these data. The flickering 
tomography could only be performed with the Ha 
line, which is more intense than the He I one. In 
Fig. 6 we show the average of all spectra used in 
flickering tomography and the average RMS spec- 
trum already corrected from instrumental, orbital 
as well as secular effects. The flickering line pro- 
file defined as above is single peaked. The RMS 
increase with the emission intensity but not as ex- 
pected from pure photon noise. It can be seen 
that the flickering and the flux line proflles present 
slightly distinct intensities in the blue side contin- 
uum region. The RMS flickering intensity relative 
to the average line intensity obtained in this work 
is lar ger than th at observed in the nova-like V442 
Oph ()Diajl200ll) . 

The Ha flickering Doppler tomogram is shown 
in Fig. 7. The filter used to produce this map 
was set to 100 km s^^ FWHM and the final tomo- 
gram resolution is estimated as 130 km s^^. The 
positions marked in the tomogram are, as in the 
flux Doppler tomograms, from top to bottom, the 
secondary center of mass, the internal Lagrange's 
point Li, the system's center of mass, and the cen- 
ter of mass of the primary. As the orbital inclina- 
tion and masses cannot be strictly determined, the 
stellar positions are plotted just as a reference. 

At least two regions of enhanced variance can 
be identified: (i) the {—Vx,+Vy) quadrant near 
the line that connect the stellar centres, and (ii) 
the {+Vx,—Vy) quadrant. Comparing this vari- 
ance tomogram with the flux tomogram (Fig. 3), 
one can see that the region of higher emissivity 
has a different shape and it is significantly shifted 
in velocity. 

The ratio of the RMS flickering tomogram and 
the flux tomogram was calculated (Fig. 5) in order 
to verify if the flickering is dependent of the line 
flux of the originating region. The high velocity 
regions of the resulting image were also excluded 
because the flux tomograms are more noisy and 



present values near zero. The resulting relative 
tomogram is then limited to intrinsic velocities be- 
low w 400 km s~^. The ratio between the tomo- 
grams does not present bright regions as do the 
flux and variance tomograms. However, the rela- 
tive flickering map is not flat but regions of lower 
intensity and a higher intensity near the center of 
the image can be seen. From this it can be con- 
cluded that the flickering can vary with the region 
of the disc where the flickering is formed, and that 
the flickering RMS activity may be independent 
of the flux intensity of its originating region. Ad- 
ditionally, the flickering in V3885 Sgr is not uni- 
formly distributed in the disc. This behavior was 
also observed in V442 Oph bv lDiad (|200ll ). 

In order to obtain flickering activity profiles, 
variance tomograms centred in the primary were 
also calculated, excluding the high noise periph- 
eral regions. Following the same method used for 
the calculation of the emission profiles, the me- 
dian along concentric rings centered on the pri- 
mary was calculated. These emissivity profiles as 
a function of velocity are shown in Fig. 8. The me- 
dian was used as the statistical criterion because 
it showed to be more robust than the average or 
mode, avoiding that the emissivity profile could 
be affected by the high noise regions. One can 
see from the upper panel in Fig. 8 that the aver- 
aged fiickering at a given velocity (or radius in a 
Keplerian disc) is roughly correlated with the disc 
emissivity for absolute velocities higher than ap- 
proximately 300 km s~^. These profiles have dif- 
ferent behavior when compared with those found 
for V442 Oph, where both the emissivity profiles 
and the ratio of the tomograms presented a de- 
pression at lower velocities. In V3885 Sgr such 
behavior is not observed. 

3.5. Flickering Simulations 

Simulations of accretion discs including a 
stochastic variability component were made, aim- 
ing to perform a better analysis of the flickering 
tomogram. The simulations were held considering 
a Keplerian accretion disc, with a line emissiv- 
ity profile given by a radial power-law. For each 
'exposure', a synthetic line profile was generated, 
considering the presence of noise with a given S/N 
ratio. Over each synthetic line profile, the contri- 
butions of isolated emission sources in the disc can 
be included to simulate the hot spot and/or the 
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boundary layer. The synthetic flux and flickering 
tomograms were generated with the resulting line 
profiles. 

The flickering activity was modeled as a series 
of flares at random positions inside a restricted re- 
gion of the disc, or at the secondary's illuminated 
face. The flickering mean frequency and ampli- 
tude are given as a function of the total line profile 
flux. These are the main simulation parameters. 
The flickering flares were included in the simula- 
tion as having instantaneous rise and exponential 
decay. 

Aiming to verify whether the flickering is pro- 
duced in the accretion disc, a tomogram of a Ke- 
plerian disc was simulated using the V3885 Sgr 
orbital period and probable masses. Comparing 
this tomogram with the observed flickering tomo- 
gram, one can verify that the position of the most 
intense flickering feature (hereafter Ci) does not 
match the accretion disc in the tomogram (Fig. 
9a). The simulations show that the absolute ve- 
locity of the flickering Ci source is not consistent 
with the velocities found in any Keplerian disk 
contained within the primary Roche lobe, more- 
over a tidal limited disk. However, the possibil- 
ity of the Ci feature being originated in a non- 
Keplerian disk cannot be excluded at all. The ab- 
solute ballistic velocities at hot spot impact region 
{vstream ~ 280 km at the tidal disk radius - see 
flg. 7) are large when compared with the Ci ab- 
solute velocity (« 130 km s^^), as expected from 
the gas accelerating in the white dwarf's potential 
well. If the actual disk radius is smaller than the 
tidal radius such a difference would be even larger. 
These two arguments suggest that Ci is not the re- 
sult of line flickering at the hot spot or outer disk 
and they are independent of the addopted orbital 
phase. 

The next step is to test the hypothesis of the Ci 
feature originated in the secondary's illuminated 
face. Synthetic flickering tomograms were com- 
puted with flickering flares from the secondary's 
inner face and these maps were compared with 
the observed flickering tomograms (Fig. 9b). We 
can verify a fairly good agreement between these 
two reconstructions. This simulation was repeated 
considering different values of stellar masses and 
orbital inclination given by the mass diagram, and 
the position agreement between the flickering fea- 
ture Ci and the inner face of the secondary star 



in the tomogram remained. 

There are at least two possibilities to explain 
the observation of line flickering at the companion 
star: first, the variability could be due to the il- 
lumination effect by flickering events produced at 
the accretion disc or, second, they could be due to 
fast chromospheric activity of the secondary. Be- 
low we will examine these two possibilities. 

To verify the plausibility of the Ci flickering 
feature being observed in the secondary due to the 
reprocessing of the flickering from the disc, some 
basic quantities are calculated. The typical recom- 
bination time-scale is given by equation 1 below, 
where iVe is the electronic density and a a is the re- 
combination coefficient. An electron density N^, « 
3x10^^ cm^'^ is need in order to observe flickering 
with a time-scale of about 1 second. 



1 



NeUA 



3 X 10^2 



10^ 

-j^yr (1) 



The recombination volume needed to produce a 
given flux F by recombination is given by equation 
2, where d is the distance to the source, ly the 
frequency of the emission and a the recombination 
coefficient. 



V 



(2) 



The flux of the Ci feature in the observed flick- 
ering tomogram was obtained by integrating the 
emission in the feature and subtracting the flux 
of the diffuse emission near it, yielding Fq-^ = 
3.6 X lO^^"* erg s^^ cm^^. In addition, a con- 
straint on the electron density is needed. A lower 
limit could be obtained from the fact that for- 
bidden lines are not observed in CVs, from this, 
the electronic density must be higher than « 10^ 
cm~^. 

Using equation 2, the lower electronic density 
limit mentioned above gives a volume oiV 10^'* 
cm'^, that is even big ger than the secondary's vol- 
ume ( Eggletoiilll983l ) so the recombination process 
could not generate the given flux in an environ- 
ment with Ne = 10^ cm~'^. 

An upper limit to the average density can be 
obtained from the clas sical atmosphere models by 
Hubenv fc La^ (|l995l ). The structure of a 10000 
K and log(g) = 4.0 atmosphere was simulated, 
resulting in an atmosphere thickness of z = 4 x 
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10"'^'' cm, measured from r^oss = 10~^ down to 
TRoss = 1- We scaled this atmosphere to T = 5000 
K using the pressure scale height, obtaining an 
atmosphere thickness z' = 1.2 x 10^ cm and a 
secondary's surface gravity log(5s) ~ 4.6 . The 
surface gravity value for the secondary in V3885 
Sgr i s consisten t with an isolated MOV or M2V 
star (|c3 I2OOOI ). Starting from the Roche lobe 
radius and the obtained atmosphere thickness, the 
average electron density is calculated (equation 2) 
as iVe « 1.4 X 10^^ cm~^. This value is an upper 
limit for the electron density, as an atmosphere 
expands as irradiated and, from equation 2, the 
electron density is inversely proportional to the 
square root of the emitting volume. Such a value 
is effectively within the optically thin region of our 
model atmosphere. 

Using high time-resolution V-band data taken 
on a single run (Bruch 2006, private communica- 
tion) we have calculated an average orbital light 
curve. The data comprise 12 orbital cycles and 
were binned into 40 phase bins using the orbital 
ephemeris proposed in this work. The average flux 
curve presents a maximum around phase 0.6 which 
is consistent with significant companion illumina- 
tion effects in the optical continuum (fig. 10). The 
possibility of the main line fiickering feature be- 
ing associated to the hot spot instead of the sec- 
ondary can not be completely ruled out from the 
data at hand. However, it would imply in a large 
(> 0.2) spectroscopic phase shift. In this case the 
broadband light curve peak could not be associ- 
ated with the secondary star. On the other hand, 
the expected absolute gas velocities at the hot spot 
are quite uncertain but possibly larger than those 
measured for the fiickering feature at 14 « 0. 

From the arguments above, we suggest that the 
Ci emission line fiickering observed in V3885 Sgr 
could be due to a radiative variable component 
produced by recombination at the secondary's illu- 
minated surface. The calculation described above 
provides only orders of magnitude, while a detailed 
model considering an irradiated atmosphere is still 
needed to ensure this conclusion. The observation 
of flickering coming directly from the inner accre- 
tion disc is not possible via flickering tomography 
because the information from higher velocities is 
diluted at the outer border of the tomograms. The 
time delay in recombination process could affect 
significantly only the flickering at high frequencies. 



As we used a integration time of about 3 minutes 
in our data, only the low frequency flickering com- 
ponent is mapped. 

The possibility of the secondary be a flare star 
was also considered. Stellar flares are observed 
through all the main sequence, with the most in- 
tens e flaring being d ominant in the hotter stars 
(see Pettersen ( 19891 ) for a review). Stellar flares 
are observed mainly in cool stars, where the flares 
have smaller amplitude and are more frequent 
than the ones of the hot stars. The mechanism 
proposed to explain the flare production involves 
the magnetic field of the star. The dynamo mech- 
anism is enhanced by the fast stellar rotation, for 
this reason it is expected that stars with higher 
rotat ion present enhanced fiare activity (| Parker 



19551 ). Cataclysmic Variables with their short or- 



bital periods are the perfect scenario for this kind 
of activity. Stellar flares were observed in S PSS 
detached binary systems ( Silvestri et al.|[2006 ) and 
also i n the magnetic CV AM Her during its low- 



state (|Kafka et al.l 120051 ) . Flickering and stellar 
fiares present some differences: first, the fiares 
appear as isolated outbursts while flickering is a 
"continuous" phenomenon and, second, the time- 
scales are dis tinct. AM Her h as an orbital pe- 
riod of 3.8 h (jYoung et al.lll98l[ ). which is shorter 
than the V3885 Sgr orbital period. Therefore, the 
secondary in AM Her has higher rotation velocity 
to feed the dynamo and presents approxim ately 



one flare at each 5.4 d ( Kafka et al. 20051 ). As 



the less intense fiares are more frequent, could the 
microfiare phenomeno n mimic the line flickering? 
Robinson et al. ( 19951 ) observed 32 flare events in 



2 h observation of the microflaring dM8e star CN 
Leo. This represents a rate of 1 flare at every 
3.7 min, which is considerably lower than the fre- 
quency of flickering events. 

The energy associated with the individual 
events in stellar flares and flickering in CVs are 
also distinct. The typical luminosity of a V3885 
Sgr line flickering event is calculated directly 
from the Doppler maps as L « 4 x 10^^ erg 
^. Using the microflare time-scale of t w 5s 



(jRobinson et al.l 11999') and the energy range of 
10^^ — 10^^ erg ( Parker! [l988f l. a luminosity range 
of 2 X 10^^ — 2 x 10^'' can be inferred for the mi- 
crofiare events. The microfiare energy range is 
therefore well below the fiickering energy. The 
fiare energy in red stars is below 10^^ erg, this. 
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combined with a typical time-scale of 10 s yields a 
luminosity below 10^^ erg s~^. This later value is 
also below the flickering flare energy found in the 
V3885 Sgr secondary star. 

As the stellar flares are distinct from the flick- 
ering both in frequency and energy, this phe- 
nomenon does not seem to be the main cause of 
the Ci emission line flickering observed in V3885 
Sgr. 

4. Discussion 

The stellar masses are constrained within an 
interval of 0.3 < Mi < 0.9 Mq and 0.25 < 
M2 < 0.55 M0. The uncertainty in Ki propa- 
gates to the derived mass ranges, so these ranges 
are just formal as well. This uncertainty does 
not include the fact that measured Ki velocity in 
CVs m ay not reflect t he w hite dwarf orbital ve- 
locity. iHartley et al. ( 2005 ) constrained the pri- 
mary mass between 0.55 and 0.8 Mq, which is 
in agreement with our mass interval. Concern- 
ing the orbital inclination, we conclude that the 
possibl e values for the sys tem are between 45° 
e 75°. ICowlev et"al] (|l977[ ) used 60° as a upper 
limit for the val ue of the orbital in c hnatio n in their 
mass diagram. iHaug fc Drechsel (|l985l) assumed 
that the orbital inclination is between 60° and 80°. 



Hartlev et al.l (|2005) pointed that the orbital incli- 
nation must be higher than 65°. The uncertainty 
on this parameter remains high. 

Regarding the morphology of Doppler recon- 
structions, an enhanced emission in the {—Vx, +Vy) 
quadrant was observed both in Ha and He I to- 
mograms. In both cases an enhanced emission 
in {+Vx,—Vy) quadrant might be interpreted as 
being originated by the impact of the gas stream 
that have not reached the disc in the {—Vx,+Vy) 
quadrant. Emission from the same region was 
searched by visua lly inspecting th e Dopp ler To- 



mogram atlas by I K ait chuck et al.l (|19941 ). The 



presence of enhanced emission in the (-1-14, —Vy) 
quadrant is seen in V1315 Aql (H/3 and He I 4471 
lines), IP Peg (H/3), LX Ser (He I 4471), DW UMa 
(H/3) and DX UMa (H/3). From all these systems, 
only IP Peg is classified as a dwarf nova, while the 
other objects are classified as nova-like CVs. The 
orbital periods of all these systems are shorter 
than 5 hours. This effect may well be due to a 
selection effect as the catalogue has only 4 objects 



with orbital period longer than 5 hours and 14 ob- 
jects with period shorter than 5 hours. High S/N 
tomography studies of other CVs are needed to 
investigate if there is any correlation between the 
orbital period or eruptive class and the presence 
of emission in the {Vx,—Vy) quadrant. 

While the H/9 a nd H7 Doppler maps presented 
by iHartlev et all (l20oi have a ring signature, 
while our Ha Doppler map is filled at low veloc- 
ities. The obtained radial emissivity indexes for 
both Ha and He I lines are similar to the one s 
obtained for V841 Oph bv lDiaz fc Ribeirol (|2003h . 
Concerning the absolute intensities, V3885 Sgr is 
approximately 4 times fainter in the Balmer emis- 
sion lines than V841 Oph. 

V3885 Sgr is the second object targeted for a 
flickering tomogr aphy s tudy, the first object was 
V442 Oph (Diazl I2OO1I ). The flickering variance 
tomogram presents two regions of higher variabil- 
ity in the {—Vx,+Vy) and {Vx,—Vy) quadrants, 
which are not in the same location as the regions of 
higher emission. In the V442 Oph variance tomo- 
gram ( Diaz 200 ih , the region of most intense emis- 
sion is in the {—Vx,—Vy) and (-1-14,-1^) quad- 
rants. The ratio between flickering and flux tomo- 
gram for V3885 Sgr and V442 Oph also presents 
a different behavior. For V442 Oph there is an 
region of low intensity near the center, while for 
V3885 Sgr the lowest intensity regions are not cen- 
tered. In both cases it is verified that the flickering 
can vary with the position within the disc. The in- 
trinsic distribution of the flickering as a function 
of the velocity also presents a distinct behavior 
for both systems. The flickering distribution along 
the V3885 Sgr disc showed to be non-homogeneous 
as observed in V442 Oph. A future study con- 
cerning methods for qualifying the significance of 
features in flickering maps is planned. 

Comparing the flickering tomogram with simu- 
lations, we have associated one of the sources of 
emission line flickering with the secondary star. 
As far as we know, this is the first detection of 
flickering at the companion star in CVs. Previ- 
ous eclipse mapping studies have located the con- 
tinuum flickering sou rces restricted to the accre- 
tion disc in UU Aqr (|Bortolettoll2"006l ) and to the 
accretion disk and the gas s tream in V2051 Oph 
(|Baptista fc Bortolettoll2002f) . Using the flickering 
tomography method we are not able to observe 
the flickering originated in the inner parts of the 
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accretion disc because the information at high ve- 
locities is spread at the outer radii of the Doppler 
map. In addition, the S/N ratio in the hne wings 
is usually low. According to our simulations, the 
brightest line flickering source in V3885 Sgr is as- 
sociated with the secondary star. One mechanism 
that may be proposed to explain this behavior is 
the flickering in the UV continuum from the accre- 
tion disc being reprocessed on the illuminated face 
of the secondary star. As the internal part of the 
accretion disc rotates faster than the disc's outer 
parts, the magnetic reconnection in the internal 
disc could be a me chanism to produce the high fre - 
quency flickering ( Geetsema fc Achterbergl 1992h . 
Flickering in the UV can be produced in the inner 
hot disc, resulting in an illumination effect in the 
optical region. The illumination of the atmosphere 
of the companion star results into both heating 
and ionization. Recombination in the optically 
thin atmosphere results in a low frac tion of the il- 



lumina tion energy by emis sion lines ( Strittmatteil 
119741) . iHartlev et al.M2005l) presented evidences of 
illumination of the secondary star in V3885 Sgr, 
as a narrow emission component superposed to the 
lines and emission from the illuminated face of the 
secondary star in the H/? and H7 Doppler maps. 
We present photometric data that supports the 
scenery of a illuminated companion star. In our 
opinion, it would be much more difficult to ar- 
gue in favor of the hot spot as responsible for the 
CI feature. The absolute velocity for the hot spot 
would be too low and it would require a large phase 
shift (> 0.2). 

5. Conclusions 

We present the analysis of an extensive spec- 
troscopic observation of V3885 Sgr. A long term 
ephemeris was determined for V3885 Sgr, with or- 
bital period Port = 0.20716071(22) d, and spectro- 
scopic conjunction, Eq = 51942.08371(46) MHJD. 
The radial velocit y semiamplitud e s and the mass 
ranges derived by iHartlev et all (|2005l ) are con- 
firmed by our measurements. 

Using ffickering tomography, it was verified that 
the flickering intensity can vary within the disc. 
From the V3885 Sgr flickering tomograms emissiv- 
ity profiles, one can see that the fiickering is most 
intense at low velocity regions. From the flickering 
tomography studies of V3885 Sgr and V442 Oph, 



one can conclude that the line flickering in CVs 
may have a non-homogeneous distribution along 
the disc. More additional flickering tomography 
studies from other systems are needed to ensure 
this conclusion. Comparing synthetic flickering to- 
mograms with observed ones, we concluded that 
the most intense line flickering source could not 
be originated in a Keplerian accretion disc defined 
by the V3885 Sgr mass interval and orbital period. 
The main flickering source may be associated with 
the secondary star, the origin of such feature may 
be the flickering in the UV continuum produced 
in the accretion disc being reprocessed at the illu- 
minated face of the companion star. 

This work is based on data obtained at LNA/CNPq 
and Cerro Tololo observatories. F.M.A.R is 
grateful from support from FAPESP fellowship 
01/07078-8. MD acknowledges the support by 
CNPq under grant #301029. We wish to thank 
Ruth Gruenwald for her carefuU reading of the 
manuscript. 
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Fig. 1. V3885 averaged spectra. The Ha and 
Hel lines are indicated. 
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Fig. 2. — Mass diagram. For a description about 
the lines see the text. The region of most probable 
masses is shown in gray. 
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Fig. 3.— Top: (a) Ha and (b) He i Doppler To- 
mograms. The central '+' marks the system's 
center of mass, above it the internal Lagrange's 
point Li and the secondary's center of mass are 
marked. Below the origin is the center of mass 
of the primary. These positions arc calculated for 
Ml = 0.57 M0, M2 = 0.42 M© and 60° orbital 
inclination. The estimated resolution for the to- 
mograms are 111 km, ,s^^ (FWHM) for Ha and 147 
km (FWHM) for He i. Bottom: Tomograms 
of (c) Ha and (d) He i emission centred at the pri- 
mary and subtracted by a radially symmetric disc 
fitting. 
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Fig. 4. — Comparison between phase maps con- 
structed from observed (a) Ha and (b) Hel 6678 
line profiles and (c) Ha and (d) Hcl phase maps 
from Une profiles reconstructed from the Doppler 
maps. Notice that the scales of the velocity axes 
are distinct for each graph panel. 
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Fig. 5. — (a) Radial cmissivity profile for Ha (tri- 
angles) and He i (circles) corrected from redden- 
ing. Simple power laws are fitted to the data, (b) 
Emissivity profiles as a fimction of the velocity of 
the emitting region, for both Ha (solid line) and 
He I (dotted). 
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Fig. 6. At the top: averaged flux (solid line) and 
flickering RMS (dotted line) spectra. In the lower 
panel the ratio between these profiles is presented. 
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Fig. 7. — (a) Flickering variance tomogram and 
(b) the ratio between the flickering and flux to- 
mograms. The '+' are, from top to bottom, the 
center of mass of the secondary, the Li Lagrange's 
point, the system's center of mass and the center 
of mass of the primary. These positions arc cal- 
culated for Ml = 0.57 Mq, Ma = 0.42 Mq and 
60° orbital inclination. The high amplitude fea- 
tures in the border of the lower panel map arc 
artifacts produced by the division by values close 
to zero (the available signal in both flux and flick- 
ering tomograms drops with increasing velocity). 
The velocity of the ballistic stream at the tidal 
disk radius is marked with a 'x'. 
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Fig. 8. — Flickering activity profiles as a function 
of the their formation velocity. The RMS curve 
is given in the top panel, while the ratio between 
RMS and flux is shown in the bottom panel. 
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Fig. 9. — (a) Synthetic accretion disk map and 
(b) synthetic tomogram of flickering coming from 
the secondary (greyscale and black contour levels) 
superposed with the observed flickering tomogram 
(red contour levels). The stellar positions shown 
in blue are the same of figure 6. 
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Fig. 10.~ V band orbital lightcurve of V3885 
Sgr. Each point is an average of approximately 
170 photometric datapoints comprising 12 orbital 
cycles. 



18 



Table 1: Journal of Observations 



observation date 


telescope 


exp. time (s) 


n. of spectra 


n. of cycles 


1999 Sep 01 


1.60-m LNA 


45 


280 


1.17 


1999 Sep 02 


1.60-m LNA 


45 


390 


1.55 


2000 Jul 07 


1.60-ni LNA 


30 


148 


0.46 


2000 Jul 08 


1.60-m LNA 


30 


72 


0.28 


2000 Jul 09 


1.60-m LNA 


30 


134 


0.54 


2001 Mar 21 


1.60-m LNA 


45 


112 


0.45 


2001 Mar 22 


1.60-m LNA 


45 


62 


0.28 


2001 Sep 10 


1.60-m LNA 


90 


128 


1.13 


2001 Sep 11 


1.60-m LNA 


90 


30 


0.31 


2001 Sep 18 


1.60-m LNA 


100 


117 


1.06 


2001 Sep 19 


1.60-m LNA 


100 


89 


0.88 


2002 Jun 16 


1.5-m CTIO 


180 


25 


0.26 


2002 Jun 17 


1.5-m CTIO 


180 


25 


0.26 


2002 Jun 18 


1.5-m CTIO 


180 


32 


0.38 


2002 Jun 19 


1.5-m CTIO 


180 


30 


0.35 


2002 Jun 20 


1.5-m CTIO 


180 


37 


0.39 


2002 Jun 21 


1.5-m CTIO 


180 


36 


0.38 


2002 Jun 22 


1.5-m CTIO 


180 


29 


0.47 


2002 Jun 23 


1.5-m CTIO 


180 


42 


0.46 
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